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I. INTRODUCTION 

The Coupled-Channels (CC) method is the most pow- 
erful tool to study nuclear reactions. However, it becomes 
extremely complicated when it is necessary to include a 
large number of channels in the calculation. In situations 
where one is only interested in a detailed description of 
a single channel, e.g. in elastic scattering, one can re- 
sort to the polarization potential approach. It consists of 
deriving a potential, to be added to the Hamiltonian of 
the elastic channel which leads to the same elastic wave 
function as that obtained by solution of the CC equa- 
tions. The natural framework to derive this potential is 
Feschbach formalism Although the exact derivation 
of the polarization potential may be as difficult as solving 
the CC problem, in some situations it is possible to find 
good approximations for it. 

Approximate polarization potentials have been derived 
by several authors, for collisions at near-barrier energies 
0, H, IU and energies well above the barrier [j| . A serious 
drawback of these potentials is that they are non-local 
and I- and .E'-dependent. Although any non-local poten- 
tial can be replaced by a trivially equivalent local one, 
the latter present poles and has an artificial dependence 
on the quantum numbers of the elastic wave function. 

For practical purposes, it is convenient to have local 
polarization potentials independent of I, in order that it 
could be used in standard computer codes. Several ap- 
proaches have been proposed to achieve this goal. In a 
recent paper, Lubian and Nunes Q tested the validity 
of the approach of Thompson et al. [7], in the case of 



breakup coupling, which is very important in collisions 
of weakly bound projectiles Q. In the present paper 
we extend this study to other procedures to derive l- 
independent potentials in the case of n Li + 12 C scatter- 
ing at near-barrier energies. For simplicity, we represent 
the continuum by an effective bound channel, restricting 
ourselves to two-channels. Although this is not an ap- 
propriate description of the continuum, it is suitable for 
the qualitative purposes of the present work since this 
effective state also leads to a long range polarization po- 
tential. 

This paper is organized as follows. In sect. 2 we give a 
brief description of the polarization potential, according 
to Feshbach's formalism. In sect. 3 we discuss differ- 
ent prescriptions suggested to derive /-independent local 
potentials without poles. In sect. 4 we apply these pre- 
scriptions to n Li + 12 C scattering and compare the re- 
sults with the corresponding ones obtained with the CC 
method. Finally, in sect. 5 we summarize the conclusions 
of this work. 



II. FESHBACH'S FORMALISM FOR 
POLARIZATION POTENTIALS 

Let us consider a collision described in terms of the 
collision degrees of freedom (projectile-target separation 
vector), r, and a set of intrinsic degrees of freedom, rep- 
resented by £. The scattering wave function 1 3K+'(r, £) 
satisfies a Schrodinger equation with the total Hamilto- 
nian 

H = H + h + V. 
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Above, H depends only on the collision degrees of 
freedom, h acts only on the intrinsic space and V couples 
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collision with intrinsic degrees of freedom. It is conve- 
nient to perform the channel expansion of the scattering 
state, in terms of the eigenstates of h, which satisfy the 
equation 



h \n) = £„ \n) . 



in the form 



(+) 



\a). 



(1) 



The intrinsic space can be divided in two complementary 
parts by the action of the projectors 



P=|0)(0|; Q = J>)(a| 



(2) 



Above, |0) stands for the ground state of h and we take 
Eo = 0. These projectors have the properties 

P 2 = P, Q 2 = Q; QP = PQ = and P + Q = 1 . (3) 

Acting with these projectors on the scattering state, one 
gets 



P 



*(+)) 



£UW) \a) = \* Q ) 



(+)\ = 



/ 



= */ 



(4) 
(5) 
(6) 



Applying these projectors on the Schrodinger equation, 
using Eq. ^ and rearranging the terms, one obtains 



[E-H] \V P )=PYQ\* Q ) 
[E-QUQ] |* Q ) = QVP|* P ) 



(7) 

(8) 



Note that PHQ = QHP = 0. 



Eq.([7J) can be used to define the polarization potential 
operator viz, 



PVQ\* Q ) = U po1 |* P ) 



(9) 



This is the basis of the calculation of Refs. [6|, [9( as they 
employ a numerical solution of the CC equations and 
make use of the above definition. To derive an explicit 
expression for the polarization potential operator, one 
first derive an expression for the projected state \^q) by 
multiplying Eq. ((SJ from the left with the Green's function 



<(+) = 1 

*W E-QMQ + ie' 



The result is 



= Gjj QYP |* P ) 



(10) 



(11) 



Inserting the above equation in Eq.([7|), we get 

[E-Meff] |*p)=0. (12) 
Above, H e ff is the effective Hamiltonian 



= PHP + PWQ Gqq QYP 



(13) 



Taking the scalar product of Eq. (|12[) with (01 , using the 
explicit form of P (Eq.©) and replacing (0| = |^ >, 
we obtain the Schrodinger equation for the elastic wave 
function in the space of the collision degree of freedom 



{K + U^ + U**) \ifo)=E \i> ). 



(14) 



Above, K is the kinetic energy operation, and the optical, 
U opt , and the polarization potentials, U po1 , are 



U opt = {0\H-K\0) 

jjpoi = ( | V q G (+) q V |o)_ 



(15) 
(16) 



The above form for the polarization potential is numeri- 
cally identical to the one defined by Eq.© 

As a consequence of the analytical structure of the 
Green operator (the presence of ie in Eq.f llOp). we can 
immediately write: 



Im {cjjU = -ttS(E- QMQ) 
while the real part is 



{ 



Re GXA = V 



QQ i ' 1 E-QUQ 



(17) 



(18) 



This last equation can be rewritten as 

which, with the help of Eq. (fT7|) gives the desired result 

Im 



Re 



{gSJ(£)} = / dz \ Q l i . , m, 



When Eqs. tfTTI) and are used for the Green function 
in Eq. ljlOp , we obtain the operator form of the dispersion 
relation. 



Re{U pol {E)} = ^-V I dz 



lm{U(z)} 
(z-E) ■ 



(20) 



One would expect the calculated polarization potential 
to satisfy the above relation. Since one usually resorts 
to several approximations to derive a local /-independent 
polarization potential, the above relation may eventually 
be broken. However, it does supply an important check 
on numerical evaluation of the polarization potential, just 
like unitarity and the S-matrix. 
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For practical purposes, it is convenient to write Eq. ([T4[) 
in the coordinate representation. While the optical 
potential is generally taken to be local and energy- 
dependent (owing to the effect of exchange non-locality, 
removed to obtain a non-dispersive local equivalent po- 
tential ) , the non-locality of Gqq leads to a non-local 
polarization potential. One obtains the equation, 



[K + U°*(t)] Mr)+I U*>° l (r,r')Mr')d 3 r' =£VoO), 

(21) 



with 



C/P°'(r,r') = (r| (0| VQ G$ Q V|0) 



(22) 



For some applications one replaces U pol (r,r') by the 
trivially equivalent local potential, 



U pol (i 



This potential has some undesirable features. Firstly, it 
has poles where ipo( r ) vanishes. Secondly, the above pro- 
cedure introduces artificial dependences on the quantum 
numbers of i/jq. These points will be discussed in further 
details in the next section. 

We should mention that our general expression for 
U pol (r,r') (Eq. ([22]) ) can account for the most general 
coupled channels situation. In particular, in the CDCC 
discussed in Ref. [6[, the breakup continuum is dis- 
cretized into an orthonormalized set of bins, which would 
then span the (J-space. The polarization potential for 
this case, after writing for Q, 



(24) 



Clearly both expressions for U pcl above are highly non- 
local by construction. In addition, the inclusion of 
continuum-continuum couplings produces further non- 
local effects @. It is easier to deal with an U po1 
having the form of Eq. (|28p . How to find an equiva- 
lent no-continuum-continuum-coupling polarization po- 
tential? To answer this question we rely on our recent 
work on the excitation of giant resonances in heavy ion 
reactions. One usually excites a given state, which itself 
is coupled to many other excited states. Using the exit 
doorway idea Q, namely the excitation of these other 
states from the ground state proceeds from the exit door- 
way (s), one is bound to attach a width to the exit door- 
ways. Labelling the doorway states by d, the Green func- 
tion takes the diagonal form 



C^(r,r') = £ 



T d {v) G«(r,r').^(r'). (29) 



J d 3 r' U pol {r,r') ^o(r')- (23) The energies 



appearing m 



td = (d\ h \d) 



1 



E-e d - (d\H\d) +ie 



(30) 



are complex. The width of any given d-state measures 
the strength of the continuum-continuum coupling. The 
above expression for U pol (r,r') with complex e^, should 
be a faithful representation the full Green function WITH 
continuum-continuum coupling. In the following we con- 
sider a much simpler two-channel case to discuss the no- 
tion of equivalent ^-independent polarization potential, 
and leave the discussion of the continuum-continuum case 
to a future publication. 



where \b) designates the b th bin, and assuming that the 
coupling is local in the r-space, becomes 



U pol (r,r')=J2 Mr) C (r, O Ty (r') 
Above, 



66' 



■>(+) 
"bb> 



(r,r') 



1 



E - QMQ + ie 



\b') 



and 



T b {v) =(0 |V(r)|6), 



(25) 



(26) 



(27) 



with an analogous expression for Tyi?' )■ If we neglect 
continuum-continuum coupling and the width of the bins, 
the Green's function becomes diagonal and Eq. (l2"5|) be- 
comes, 



U po1 {y,y') = ^2 Mi 



<(+) 

"6 



(t,t>) F b {r> ). (28) 



One frequently perform angular momentum projec- 
tions in the Schrodinger equation. For the simple case 
of a scalar Hamiltonian with scalar coupling, one makes 
the expansions 



(31a) 



U pol (r,r') = ± 7 J2Y lm (r) Uf°\r,r') Y lm (v')- (31b) 



The angular- momentum projected version of Eq. (|23[) is 
the local equivalent polarization potential, 



U[ ol (r) 



ui(kr) 



dr' Uf°\r,r')ui{kr'). (32) 



This potential should be included in the Hamiltonian for 
optical model calculations of the elastic radial wave func- 
tion, ui (r) . It this way, one needs a different polarization 
potential for each partial wave. 

The polarization potential of Eq. (|3"2")l is not very use- 
ful, since it requires the knowledge of the exact radial 
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wave function in the elastic channel, ui(kr). Determining 
ui(kr) is as hard as solving the original CC equations. 
A rather widely used approximation for the /-dependent 
local polarization potential consists of replacing ui(kr) 
by the optical radial wave function, wi(kr), which is the 
solution of the partial-wave projected Schrodinger equa- 
tion with the polarization potential switched off. This 
potential, denoted by Ui, is 



Ui(r) 



wi{kr) 



dr' U[ ol {r,r') wi{kr'). 



(33) 



A further inconvenience of trivially local equivalent po- 
tentials is that they have poles wherever the radial wave 
function vanishes. In what follows we discuss three dif- 
ferent prescriptions for obtaining /- independent polar- 
ization potentials free of poles. 



The modified Thompson prescription (MT) above guar- 
antees the contribution of the /— values around the graz- 
ing one, regardless to the behavior of |u;(r)| 2 . One phys- 
ical motivation for choosing dTi/dl instead of 7] is that 
the DWBA amplitudes of non-elastic processes in the adi- 
abatic limit do behave as the /— derivative of the elas- 
tic S'-matrix elements, which enters in the definition of 



C. The semiclassical prescription 

This prescription relies on the semiclassical idea that 
the orbital angular momentum, if treated classically, 
should be related to r and the energy E through the 
definition of the classical turning point r t , namely 



III. LOCAL, X- INDEPENDENT 
POLARIZATION POTENTIALS 

A. The prescription of Thompson et al. 

Thompson et al. Q proposed the following definition 
of an /-independent version of the polarization potential, 



£ ; (2/ + l)r, Mfcr)| 2 u?° l (r) 
(2/ + 1)7} | Ui (At)P 



U T {r) 



(34) 



where TJ is the transmission coefficient in the elastic 
channel for the I th partial-wave. The above definition 
of U T guarantees that no poles remain in the polariza- 
tion potential, which arise from the presence of ui (r) in 
Eq. (pM)) . Furthermore, owing to the presence of TJ in 
the /-sum, only values of / where Tj is close to unity will 
contribute. Of course the probability density |u/(r)| 2 is 
small for small values of / due to absorption. Thus, in 
the prescription of Thompson et al., Eq. (|34[) should con- 
tain contributions of / in the vicinity of the grazing one. 
We doubt that this is guaranteed always since there is 
interference effects in |w;(r)| 2 which may end up allow- 
ing the contribution of small / as well. This prescription 
has been recently used in ref. [6j in the context of the 
Continuum Discretized Coupled Channels calculation of 
break up and elastic scattering of 8 B. 



B. A modified form of the prescription of 
Thompson et al. 

We now introduce a slightly modified version of the 
above discussed prescription, by substituting T/ by its 
derivative with respect to /, namely, 



U MT {r) 



(2/ + f) (gyg) \ui(kr)\ 2 Uf ol {r) 
£, (2/ + 1) (dTi/dl) |u,(fcr)| 5 



2/i n 



l{l + l) + U{r t ) = E, 



(36) 



where U is the real part of the optical potential contain- 
ing the nuclear and the Coulomb pieces, ft is clear that, 
for a given collision energy, r t = r t (l) is a function of 
/. The prescription consists of identifying r = r t (Z) in 
order to build the /-independent polarization potential, 



U^(r) = Ur(r t (l)). 



(37) 



Using Eq. ([3"T]) . one gets the potential U sc at a discrete 
set of r- values (one for each partial- wave) and interpolat- 
ing between these points one obtains a continuous func- 
tion. 

All three /-independent potentials discussed above con- 
tain a further energy dependence, besides the one that 
arises from the green function in the polarization poten- 
tial. This extra energy dependence is non-dispersive and 
thus could render the applicability of the dispersion re- 
lation questionable. We shall verify this point in the 
following section. 

In all of the above prescriptions the starting point is the 
trivially-equivalent local potential obtained from angular 
momentum projected versions of Eqs. (|22p and (123|) . or 
the solution of the CC equations used in Eq.©. This 
will be shown in detail in the next section. 



IV. 



APPLICATION IN THE n Li+ 12 C 
SCATTERING 



(35) 



We have chosen the system 11 Li + 12 C to test the dif- 
ferent local approximations for the polarization poten- 
tial discussed in the previous sections. We have decided 
not to consider the coupling to the continuum explic- 
itly, but rather representing the continuum by a single 
bound effective channel. It is clear that this is not an 
appropriate representation for the continuum. One ob- 
vious shortcoming of this model is that it does not con- 
tain continuum-continuum couplings, which is known to 
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play an important role in nuclear reactions with weakly 
bound projectiles Q. The justification for considering 
this schematic application is that our purpose is to com- 
pare different approximations for the polarization poten- 
tial, rather than to perform quantitative calculations of 
elastic or fusion cross sections. 

The two CC equations used in our calculation can be 
easily readout from Eqs. ([7]) and (J5J. We denote the 
two intrinsic states by |0) (ground state) and |1) (ex- 
cited state). For simplicity, the numerical calculations 
are performed in the sudden limit. In this limit, the 
excitation energy of the continuum states are neglected 
(eo = £1 = 0). The projectors are P = |0) (0| and 
Q = |1) (1| and , after taking the optical potentials in 
the two channels to be the same, Eqs. ((TJ) and (J8]) be- 
come 

[E-K- U opt (r)] # (r) = ^(r) *i(r) (38) 
[E-K- U opt (r)] * x (r) = F(r) * (r). (39) 

Above, !F(t) is the complex and symmetric form factor 



T{v) = / d 3 x ^o(x) V(r,x) ^(x), (40) 

evaluated with the coupling interaction 

V(r,x) = U hT (v hT ) + U hT (v hT ) - U° pt (r). 

Since the intrinsic states are orthogonal and U° pt (r) does 
not depend on x, Eq. (|4"0")) reduces to 

J d 3 x<p (x) [Uf lT ( r hT) + U hT {r hT )] <px(x). 

For the present application, we assume that 11 Li breaks 
up into two fragments. The first, f\, corresponds to a 
neutron pair, which we treat as a single particle (di- 
neutron). The other, /2, is the 9 Li— core. In this way, 
the coordinates appearing in Eqs. (j40|) and (jITj) are 



r /lT = r + 7 1 x, 7! 
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r /2T = r + 7 2 x, 7 2 = - — . 

and the projectile-target separation vector, r. 
[//^(r/jr), t r / 2 T(r/ 2 r) and U opt {r) are the corre- 
sponding interactions. Of course T(r) can be evaluated 
exactly numerically once these potentials are given and 
the single particle wave functions of the halo neutron 
in the ground and in the excited state, <po( x ) = ( x |0) 
and </?i( x ) = ( x |1) ) are used to evaluate the integral in 
Eq. (|40"|) . For our purposes we parametrize T(r) as 



T{y) = exp 



r 
7i a 



with 



y/2m- 2 B 



(42) 



(43) 



Above, /ii-2 is the reduced mass of the fragments inside 
the projectile and B is the breakup threshold. For 11 Li, 
/ii_2 = 18 mo/11, B = 0.376 MeV, and one gets a — 
5.83 fm. For the purpose of simplicity, we take T§ to be 
real. We also ignore the Coulomb coupling altogether. A 
qualitative justification for the approximation of Eq. ([42|) 
is given in the appendix. 



The optical potentials can be written as 
U opt (r) = U c (r) + U N (r), 

where 

Z P Z T e 2 

U c (r) = ' 



2R C 

Z P Z T e 2 



o _ _r_ 



and 

U N (r) 



To 



for r < R c 
for r > R c 

-Wo 



1 + exp [{r — R r ) /a r ] 1 + exp [(r — Ri) /ai\ 

Above, R c = r oc (a 1 / 3 + A^ /3 ) , R r 

r 0r (A P /3 + aJ/ 3 ) and R, = r 0l (a 1 ^ + A^ /3 ) . In the 

above equations, we use typical values for the parame- 
ters: 



To = 3.0 McV; r oc = 1.4 fm 

V = 60 MeV; W = 60 MeV 

?*0r = 1-25 fm ?*oi = 1.00 fm 

a r = 0.60 fm a* = 0.60 fm. 



(44) 



The reduction of the above equations to get the po- 
larization potential proceeds as in the previous section 
and we get for the Schrodinger equation for ipo(v) the 
following 



[E - K -U opt (r)-U pol (rj\ Vo(r) = 0, (45) 



with 



U pol {r,r')=T{r) Gl +) (r,r') T{r'). (46) 



In order to perform numerical calculations it is con- 
venient to carry out the usual partial- wave expansions. 
Since we are using scalar form factors, the partial-wave 
projected polarization potential is given by 

U pol (r,r / ) = T(r) Hr')- (47) 

The partial-wave projected optical Green function in 

channel- 1 can be written [l2j 

= -J^ e- i5 > Wl (k ir< ) ^ (+) (fcir>). (48) 
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Above, fj, is the reduced mass of the projectile-target sys- 
tem, k 1 = y / 2/ZBi/fi ) r< (r>) is the smaller (larger) of 
the radial separations r and r' , wi(k±r) is the regular so- 
lution of the optical Schrodinger equation (partial-wave 
projected Eq. (j39|) setting T(v) = 0) with the asymptotic 
form 



wi{kir — > oo) 



H\ \k 1 r)-Si H\ >(k ir 



(-)/ 



(49) 



Above, H^\H^ +) ) 



IS 



and <5; is the nuclear phase-shift 
the Coulomb wave function with ingoing (outgoing) 
boundary condition and Si = exp (2« Si) is the nuclear 

S-matrix at the I th partial-wave. In Eq.gH]), 
is the solution of the same optical Schrodinger equation, 
but with a different asymptotic behavior. At large sepa- 
rations, it is the outgoing wave 



n{ +) (k ir oo) 



H[ + \k iroo ). 



To determine at finite projectile-target separations, 
the radial equation must be numerically integrated in- 
wards, starting from a large r— value where the above 
asymptotic form is valid. This can be easily achieved 
using a negative mesh step in any conventional code for 
numerical integration. 

The dispersion relation is clearly satisfied by Eq. (|4T|) . 
Is it satisfied by the I — dependent potential Ui (inserting 
Eq. ljlT)) in Eg. ([52)0 ? We remind that the expression used 
to get Ui is Eq. ([3"3")) . One would expect the dispersion 
relation to be satisfied by Ui only if ui(kr) is real, which 
is certainly not the case. However, a remnant of the 
dispersion relation should still be seen in Ui as has been 
shown over and over again in the study of the Threshold 
Anomaly Not withstanding the above reservations 
we will give below an account of our calculation of the 
equivalent /-independent potential. 

In figure 1 we show polarization potentials at the col- 
lision energies (a) E = 2 MeV, (b) E = 4 MeV and 
(c) E — 6 MeV. We present results obtained using 
the Thompson (T) prescription (Eq.fJM))), the modified 
Thompson one (MT) ( Eg. ([55)1 ) and the semiclassical pre- 
scription (SC) (Eq.([37])). In the cases of the T and the 
MT potentials, we evaluate the radial wave functions and 
the transmission coefficients solving the CC equations 
and then carry out the I — averages of Eqs. ([3"4")l and ([55)) . 
Also shown are the approximate Z-dependent potential Ui 
fEq. ([3"3")0 for I = and / = 10 and the optical potential 
employed in our two coupled channels. The singularities 
of Ui, arising from the nodes of the radial wave functions, 
are treated by the technique developed in Ref. [liij ]. We 
see clearly that whereas the Ui=q is quite oscillatory, the 
T, MT and SC ones seem to behave smoothly in so far 
as the real part of the potential is concerned. Oscilla- 
tions in Ui occur near the nodes of wi. Since this wave 
function is complex and their real and imaginary parts 
do not vanish simultaneously, the polarization potential 
remains finite. 
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FIG. 1: Real and imaginary parts of the polarization poten- 
tials (in MeV) calculated using the different recipes discussed 
in the text. T: Thompson, MT: modified Thompson, SC: 
Semiclassical, and the I — dependent potentials for I = and 

I = 10. The coupling stength is taken to be 3 MeV, a) 
E = 2 MeV, b) E = 4 MeV and c) E = 6 MeV, for the system 

II Li + 12 C. See text for details. 



In figure 2 we show the results for the elastic scatter- 
ing angular distributions at E = 2, 4, 6 MeV. The coupled 
channels results are shown as the full circles. None of the 
/-independent polarization potentials seems to work very 
well. However, the T potential is better then the others. 
It is close to the CC results, except for the collision energy 
of 6 MeV. In this case the angular distribution obtained 
with this potential oscillates out of phase with respect 
to the correspondint CC results. On the other hand, the 
approximate Z-dependent potential Ui reproduces accu- 
rately the CC results at the three collision energies. 

In figure 3 we show our results for the fusion cross 
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FIG. 2: The ratio of the elastic angular distribution to 
Rutherford calculated with CC (coupled channels), OM (op- 
tical model - no coupling) and with the / — dependent po- 
larization potential, Ui. The other symbols are the same as 
figure 1. See text for details. 



section excitation function, obtained by the relation 



a F (E) 



k 
E 



(Vo|-Im{[/ P t }|Vo 



The wave function tJjq was calculated solving the 
Schrodinger equation for the elastic channel, including 
the optical and each of the above discussed polarization 
potentials. Here, both Ui and the T-potential seem to 
work rather well over the energy range considered. The 
SC overshoots in the barrier region but otherwise it re- 
produces the CC calculation. The MT seems to be quite 
off, however on the average it works well too. 

Finally, to be sure about the consistency of our calcu- 
lation we have checked the dispersion relation. We have 
fixed r = 12 fm and plotted in figure 4 the resulting be- 
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FIG. 3: The fusion excitation function for the same cases as 
in figure 2. a — <™ + *— 
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FIG. 4: The energy dependence of the real and imaginay parts 
of the polarization potential (in MeV), indicating the trend 
dicatated by the dispersion relation. See text for details. 



havior of the real and imaginary parts of the polarization 
potential vs. E. It is clear that the Ui for I — shows the 
general trend of what one would expect from the disper- 
sion relation: maxima in the real part accompanied by 
sharp variations in the imaginary part. Similar behavior 
is found for the T potential. Note the unphysical sharp 
oscillations of the SC potentials. They arise from the fact 
that this potential can only be evaluated over a sparse 
mesh of r-values,correponding to the turning points at 
each angular momentum. Thus the T and U\ potentials 
exhibit the threshold anomaly In particular, the 

T potential shows a maximum in the imaginary part at 
E = 3.9 MeV, whereas the real part shows a maximum 
at E — 3.5 MeV. This behavior of the T-potential is in 
line with the breakup threshold anomaly of Ref. [l5[ • 



V. CONCLUSIONS 

In this work we have investigated I — independent po- 
larization potentials in a schematic two-channel model, in 
which the range of the polarization potentials is equiv- 
alent to that associated with the breakup channel. We 
used several prescriptions to derive I — independent po- 
larization potential and found that none could reproduce 
satisfactorily the results of coupled channel calculations. 
This conclusion is consistent with that of realistic CDCC 
calculations ((J using the Thompson prescription. On 
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the other hand, we have shown that an approximate 
I — dependent potential, obtained using the unperturbed 
Green's function in Feshbach theory, gave reasonable de- 
scriptions of elastic angular distribution at energies above 
the barrier. 



where t stands for the cosine of the angle between the x 
and the r vectors. The main contributions to this integral 
comes from 7X anti-parallel to r. We than replace 

6 (|r + 7x| — Rr) — » (r — \j\ x - R T ) 



VI. APPENDIX 

We give below, a qualitative justification of the 
adopted parametrization of F{r). Since the diffusivity of 
the potentials is very small as compared to a, R\ and R2, 
we set a r ~ 0. In this way, the Woods-Saxon potentials 
take the forms of the step functions, 



U fT {rf T ) a -V 6 (|r + 7 x| - R T ) ■ 



(50) 



Above, / = fi or /2 and 7 = 71 or 72. 

At large projectile-target separations, r > Rp + Rt, 
only x > Rp contributes to the integral of Eq . (|4"Tj) . Since 
we are neglecting angular momenta, in this range the 
ground state wave function depends only on the radial 
and has an exponential form. That is, 



tp {x) oc exp^-^j, 



with a given by Eq. (|43|) . Since we are adopting the sud- 
den approximation, the wave number of the state ipi is 
supposed to be vanishingly small. Therefore, this wave 
function is constant within the integral. The contribution 
from the potential to the form factor takes the form 

T{r) cx J dx exp J 6 (|r + 7x| - R T ) dt, 

(51) 



and assume that the integration over t does not depend 
strongly on r, leading only to a renormalization of the 
strength of the form factor. The integral of Eg. ([51} then 
becomes 



T{r) cx 



(r-R T )/\l\ 



x\ 


r 


— cx exp 




a.) 



Of course, the form factor will be dominated by the con- 
tribution from the di-neutron, which has a longer range 
(I71I > I T2 1 ) - Denoting by the constant of proportion- 
ality, we can write 



T(r) ~ T exp 

7i" 



which is adopted in the present calculation. 
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